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ABSTRACT 



Context. Maser emission from the H2O molecule probes the warm, inner circumstellar envelopes of oxygen-rich red giant and supergiant stars. 
Multi-maser transition studies can be used to put constraints on the density and temperature of the emission regions. 

Aims. A number of known H2O maser lines were observed toward the long period variables R Leo and W Hya and the red supergiant VY CMa. 
A search for a new, not yet detected line near 475 GHz was conducted toward these stars. 

Methods. The Atacama Pathfinder Experiment telescope was used for a multi-transition observational study of submillimeter H2O lines. 
Results. The 533 - 440 transition near 475 GHz was clearly detected toward VY CMa and W Hya. Many other H2O lines were detected toward 
all three target stars. Relative line intensity ratios and velocity widths were found to vary significantly from star to star 

Conclusions. Maser action is observed in all but one line for which it was theoretically predicted. In contrast, one of the strongest maser lines, 
in R Leo by far the strongest, the 437 GHz 753 - 660 transition, is not predicted to be inverted. Some other qualitative predictions of the model 
calculations are at variance with our observations. Plausible reasons for this are discussed. Based on our findings for W Hya and VY CMa, we 
find evidence that the H2O masers in the AGB star W Hya arise from the regular circumstellar outflow, while shock excitation in a high velocity 
flow seems to be required to excite masers far from the red supergiant VY CMa. 

Key words. Stars: AGB and post-AGB - Stars: individual: VY CMa, R Leo, W Hya - supergiants - circumstellar matter 
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1. Introduction 

Water (H2O) is one of the most abundant molecules in the 
atmospheres and envelopes of oxygen-rich red giant and red 
supergiant stars (see, e.g., Tsuji 1964; Ohnaka 2004). Often 
very intense H2O maser emission in the 22.2 GHz /a-„a:J1] 
= 616 - 523 transition has been found toward many hundreds 
of low- and high-mass star-forming regions and circumstellar 
envelopes of Mira and semi-regular long-period variable stars 
(LPV s), as well towar d a few luminous red supergiants (see, 



e.g., Valdettaro et al. Ibool ). Since the 616 rotational energy 



level is high (643 K) above the ground-state and its transition 
probability is rather small (because of its low frequency) this 
line is readily observable even from sea level. 

In addition to the 22.2 GHz line, maser emission from 
a number of other ortho and para H2O (sub)millimeter- 



Send offprint requests to: K. M. Menten 

' Energy levels of the asymmetric rotor molecule H2O are denoted 
JKa.K^ where J is the total rotational quantum number and K,, and 
are its projections on the orthogonal symmetry axes of the limiting 
prolate and oblate symmetric rotor, respectively. For ortho transitions 
Kg + Kc is odd and for para transitions Ka + is even. 



wavelength transitions from within the vibrational ground state 
has been detected toward sta r-forming regi o ns a nd red g iant 



and s t tpergiant stars ( Mente n et al.lll990allbt ICern icharo et al 



'1990"; iMelnick et al]|l993t I Yates et alJll995t iLiliestrom et al 
i2002i) . All of these lines, except for the 183 GHz 3i3 -22o 



transition ( Cernicharo et alJ [l990). were also observed in the 
present study and are listed in Table Because water vapor 
in the Earth's atmosphere absorbs heavily in many of these 
(and other) H2O lines, they are only observable from high, 
dry sites and several only during the small percentage of time 
when the atmosphere's H2O vapor content is minimal over the 
Atacama Pathfinder Experiment's 5 100m high site on the Llano 
de Chajnantor; see Fig. [1] Note that, in particular, the newly 
discovered 475 GHz line must be amongst the submillimeter 
lines with the poorest transmission detected from the ground 
so far. 

Toward the hot, innermost envelopes of red giants 
and supergiants rotational lines from within the lowest 
vibrationally exc ited state, the bending mode (V2-I), have 
also been found jMenten & Melnicklll989l; iMenten & YoungI 



- Frequencies and energy values are taken from the JPL catalog: 
|http://spec.jpl.nasa.gov/ftp/pub/catalog/cO 18003. cat 
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Fig. 1. Calculated atmospheric zenith transmission for the site 
of the APEX telescope for three values of the precipitable 
water vapor column: 0.2 mm (upper curve), 0.4 mm (middle 
curve) and 0.8 mm (lower curve). The upper panel shows the 
frequency range around H2O lines that have been observed 
with APEX in the 870 window, while the lower panel 
shows the frequency range around lines in the 650 /im window. 
Positions of detected H2O lines are marked as vertical lines, 
H2O lines that were searched, but not detected as dashed 
lines. Observations of the 293, 297, an d 336 GHz lin e s from 
the V2=l state have been reported by iMenten et al. (l2006h . 



observations of all the other lines in the present paper. Note 
that, both, FLASH and APEX 2a are double sideband receivers 
with intermediate frequency bands centered at 3 and 6 GHz, 
respectively. Signal and image sideband are thus separated by 
6 and 12 GHz, respectively, and frequently have significantly 
different transmission. 



wide range of physical conditions. These authors do not predict 
maser emission for the 437 GHz 753 - 650 transition first 



detected by iMelnick et alJ (Il993h in evolved stars and also 
observed by us, but do predict maser action in the extremely 
high excitation 17443 - 167,10 line near 355 GHz line, which 
we searched for, but did not detect. 

Using the Atacama Pathfinder Experiment (APEXH) 12m 
telescope we conducted a multi-transition submillimeter H2O 
maser study, observing the lines listed in Table [1] Our dataset 
is complemented by spectra of the 22.2 GHz H2O line taken 
with the Effelsberg 100m telescope close in time to our APEX 
observing run. Our targets were two asymptotic giant branch 
(AGB) objects, the Mira star R Leo and the semiregular (SRa) 
variable W Hya, and the red supergiant (RSG) VY CMa. 
The f ormer two are typical, nearby, LPVs with "reprocessed" 
(Knap p et al.ll2003h Hipparcos parallax distances, D, of 82;!;^' 
pc and 78 +6 pc, respectively. Their mass loss ratefl M, derived 
from CO observations are several times 10"^ Mq yr ' ( YoungI 



19951), although for both objects si gnificantly different values 



have been derived (for R Leo, see iMenten & Melnickl 11991 
and ESlfor W Hya). 



a mass-loss rate of 1-2 10 '*Mo yr ' (IDanchi et al 



With a luminosity, L, of 2 IQr' Lq ( Sopka et al 
^0 yr 



19851) and 



1 1994 the 



peculiar red supergiant VYCMa is a remarkable object by 
any standard and the most prolific stellar (OH, H2O, and 
SiO) maser so urce known. Its di stance is usually assumed 
to be 1.5 kpc dLada & Reidlll978h . However, a recent VLBI 
trigonometric parallax of VY CMa's SiO masers yields a 
distance a; 20% closer (Reid & Menten, in preparation). In the 
following we shall assume that the distance is 1.1 kpc. 

Our observations are summarized in ^ In ^we describe 
our results. In particular, we report the detection of a new 
maser line, the 533 - 440 transition near 475 GHz. We 
discuss the observed spectra, their similarities, but also very 
pronounced differences in the relative luminosities and the 
spectral appearances of the lines for each object and also 
among objects. In ^our observational results are compared 
with predictions of calculations that aim at modeling H2O 
maser emission. We present a summary in ^ 



2. Observations and data reduction 



19951; iLiliestrom et alj l2002h : see iHunter et al.l (l2007h for a 
compilation and interferometric measurements. Given their 
excitation conditions, these lin es must arise from w ithin a few 
stellar radii of the photosphere jMenten et al ] l2006l) . 

hereafter NM91) performed 



Neufeld & Melnickl (11991 



statistical equilibrium calculations to analyze the excitation 
of H2O masers in the vibrational ground state and outlined 
parameter regions (of, essentially, temperature and density) for 
which the energy levels of certain transitions become inverted. 
They predicted inversion in all the lines listed in Table[Tl except 
for the 437 and 443 GHz lines. Maser action has indeed been 
observed in all but one of the predicted lines toward star- 
forming regions and/or red giant stars. (Maser emission in 
the 475 GHz 533 - 440 is for the first time reported here.) In 
particular, NM91 find the very widespread 22.2 GHz 616 - 523 
transition to be the most strongly inverted one over a very 



Our observations were made between 2006 June 7 and 19 under 
generally excellent weather cond itions with the 12ri i Atacama 
Pathfinder Experiment telescope jGiisten et al ] l2006l) . 

All lines considered here are listed in Table[T] The 434, 
437, 439, 443, 471, and 475 GHz lines (see Table [B were 
observed with the First Light APEX Submillimeter Heterodyne 
instrument (FLASH, Heyminck et al. 2006) while the APEX 
2a facility receiver (iRisacher et alJl2006 ) was used to observe 
the 321, 325, and 355 GHz lines as well as the 7 = 3- 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 

* In the following all quoted mass-loss rates and luminosities are 
values taken from the quoted references, but corrected for the revised 
distances quoted in the text. 
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Fig. 2. Energy level diagram of para H2O (left) and ortho H2O (right). Upper and lower levels of known maser lines appear in 
bold, are connected by arrows and have the rounded transition frequency (in GHz) indicated. Data for all these lines (except for 
the 183 GHz line) are reported in the present paper. 



2 line of CO. Calibration was obtained using the chopper 
wheel technique under consideration of the very different 
atmospheric opacities in the signal and image sidebands of 
the employed double sideband receivers (see Fig. [TJ. The 
radiation was analyzed with the MPIfR Fast Fourier Transform 
spectrometer, which provides 16384 frequency ch annels over 
the 1 GHz intermediate frequency bandwidth (IKlein et al 



I2OO6 I). To increase the signal to noise ratio, the spectra were 
smoothed to effective velocity resolutions appropriate for the 
measured linewidths, typically ~ 0.5-1 km s To check the 
telescope pointing, the receiver was tuned to the 437 GHz 
H2O line, which showed strong emission in all of our three 
sources, and five point crosses centered on the stellar position 
with half beamwidth offsets in elevation and azimuth were 
measured. Pointing corrections were derived from the latter 
measurements. The pointing was found to be accurate to within 
«i 3", acceptable given the FWHM beam size, 0b, which is 20" 
FWHM at 321 GHz and 13" at 475 GHz. 

Additional data had been taken earlier toward VY CMa for 
the 355 GHz line (see Table[Tll for which maser action had been 
predicted by excitation modeling (see l4.4l l. These observations 
were made in 2005 July/August. 

In Table |2] we present our line intensities in a flux density 
scale (i.e., in Janskies) assum ing the aper t ure ef ficiencies 
observationally determined by iGiisten et al. (2006) for the 
respective frequency ranges. 

Observations of the 22.2 GHz 616 - 523 transition were 
made with the EfFelsberg 100m telescope on 2006 June 27, i.e.. 



ca. 1-2 weeks after the submillimeter observations. The line 
was detected with the facility high electron mobility transistor 
receiver and autocorrelator backend. The data were corrected 
for atmospheric opacity and variations of the telescope's gain 
curve with elevation. 

Fig. |2] shows our observed H2O lines on energy level 
diagrams of para and ortho water 



3. Results 

3.1. General properties of the observed emission - 
luminosities 

For W Hya, the observed isotropic photon luminosities of the 
detected H2O maser lines are all within one order of magnitude, 
ranging from 3 10"^' to 3 10^^ s"' (Tableland Fig.©. This is 
also true for VY CMa, but here the much higher luminosities 
range from 6 lO'*'* to 5 10"*^ s"^ The situation is different with 
R Leo. Here the 321 and 439 and the marginally detected 471 
GHz line have comparable luminosities of 8 lO'*', 4 lO'*' s"', 
and 4 10'*' s respectively. However, the 22.2 GHz line is 
anomalously weak (3 10"*" s"') and the 437 GHz line extremely 
bright (3 lO"*^^ s-i). 

We emphasize in particular that for W Hya and VY CMa 
the luminosity of the 22.2 GHz lines is comparable to that of 
other lines and that for R Leo it is much lower This is important 
for the discussion of maser excitation (see 94.1b . 
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Fig. 3. Top to second panel from bottom: H2O spectra observed toward VY CMa, W Hya, and R Leo (left, middle, and right 
column, respectively). The different H2O hnes are labeled. In the VY CMa panels also the approximate rest frequencies and 
energies above ground are listed. To bring out weak features and better show the full extent of its emission, the spectrum of the 
22.2 GHz line is also shown magnified by a factor of 100. The bottom panels shows the 345.8 GHz 7 = 3-2 line of CO. Note 
that some of the spectra have very poor signal-to-noise ratios, which may lead to "broadening" of the feature; cf. the 439 and 475 
GHz lines in W Hya. Similarly, it is difficult to tell how much of broad emission seen in some of the lines from V Y CMa is due 
to blended naiTow features. Also, the "bump" at v > 45 km s ' in the 437 GHz spectrum of W Hya is almost certainly a baseline 
artefact. These facts should be kept in mind in the discussion in ^3.11 



3.2. Spectral appearance - clues to the emission 
regions 

As can be seen from Table [T] and Fig. [21 the lines observed 
by us arise from energy levels that span a very wide range of 
energies above the ground state, with the lower level energy, 
E(/k, ranging from 454 to 1846 K. The highest excitation 
lines (437 GHz/1505 K and 321 GHz/1846 K) cover a smaller 
velocity extent than lower excitation lines. One might expect 
this if they are formed in a hot region closer to the star with an 
outflow velocity that is smaller than further out from where 



the lower excitation lines are arising. However, there is no 
clear monotonous trend from larger to smaller linewidth with 
increasing energy above ground. 

All the H2O lines observed toward R Leo and W Hya (Fig. 
|3]l show a single feature or a blend of features centered on or 
within a few km s ' of their stellar velocities determined by the 
CO line centroids. The features are asymmetric, except for the 
437 GHz line toward R Leo, which is nearly Gaussian with a 
FWHM of 0.78 km s ^ In contrast, for VY CMa essentiaUy 
every line has a diff'erent shape. 
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Table 1. Ground-state water maser lines observed with APEX. Table 2. Results of APEX observations 
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Frequency 


E[/k 
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J" = 


(MHz) 


(K) 


616 


- 523 


22235.08 


642.5 


IO29 


-936 


321225.64 


1845.9 


5|5 


- 422 


325152.92 


454.4 


174.13 


- 167,10 


354808.9 


5764.3 




- 660 


437346.67 


1503.7 


643 


- 550 


439150.81 


1067.7 


752 


- 661 


443018.30 


1503.7 


642 


-551 


470888.95 


1041.8 


533 


-440 


474689.13 


702.3 



Line 



Columns are (from left to right) quantum numbers of upper and 
lower state, frequency and energy above ground of lower state in 
Kelvins; k is the Boltzmann constant. Frequency values, taken from 
the JPL catalog, have formal uncertainties of order 50 kHz. More 
accurate va l ues fr om a fit to the H2O spectrum have been presented by 
IChen et alj fcood) . The difference between their values and the ones 
used by us is typically of order 20 kHz, corresponding to less than 
0.02 km s"', which is smaller than the uncertainties in our velocity 
determinations. 



For both R Leo and W Hya, the FWZP width of the widest 
H2O line is only x 20% of the CO line's width. In contrast, 
for VY CMa, some lines cover « 50 km s"', which is about 
half the velocity range of the CO line. This indicates that for 
all three stars the H2O maser emission arises from a part of the 
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Fig. 4. Ratios of the integrated submillimeter H2O maser line 
flux densities (marked by their frequency rounded to GHz) 
to the integrated line flux density of the 22.2 GHz line for 
our three program stars, R Leo (filled circles), W Hya (open 
circles), and VY CMa (star symbols). For upper limits, 3 times 
the rms levels given in Table[2]were used. Note that Ix for a line 
is proportional of its photon luminosity. 
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Columns are (from left to right) line ID (i.e. frequency in GHz; 
see Table [T), integrated flux density, velocity range covered by line 
(FWZP), peak flux density of strongest emission, velocity of strongest 
emission and isotropic photon luminosity. Numbers in parentheses 
following integrated flux densities are formal errors derived from 
Gaussian fitting and/or noise level analysis. The overall calibration 
has an estimated uncertainty of 20%. For non-detections, the quoted 
flux density (in parentheses) is three times the rms noise value of the 
spectrum smoothed to x 0.5 km s"' resolution and the quoted value of 
the integrated flux density has been calculated assuming a width equal 
to that of the H2O line found in the respective source that covers the 
widest velocity range. To calculate isotropic luminosities, distances 
of 82, 78, and 1100 pc have been assumed for R Leo, W Hya, and 
VY CMa, respectively. For upper limits of the latter three times the 
rms noise in the integrated flux density quoted in column 2 has been 
assumed. 



envelope in which the outflow has not yet reached its terminal 
velocity (i.e., half the FWZP linewidth of the CO line shown 
in the bottom panels of Fig. O. This is plausible, given the 
small size emission region of the 22.2 GHz emission, which 
for W H ya was found to be a ring with a diameter of 0."3 
(24 AU, 



Reid &Mentenlll99 0) and for VY CMa has a size 



of 0."7 X 0."4 t hat translates into a m uch larger linear size of 
770 X 440 AU jRichards etalJll998h . The latter is still much 
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sma ller than size of the CO emission region of several thousand 
AU dMuUer et al.ll2007h . 



In contrast, a spectrum of the thermally excited 557 GHz 
lio - lot ortho ground-state H2O line observed toward W Hya 
with the Odin satellite has a wid th (and centroid veloci ty) that 
is consistent with the CO values dJusttanont et al ] |2005h . 

For VY CMa the situation is more complicated. Toward 
this source, various (but not all) thermal and maser lines 
from SiO show velocities equal or even exceeding the values 
seen for CO, indicating high velocit y motions that do not 
fit into a continuo us outflow picture (ICernicharo et alJ 1 1 997 : 
Herpin et alJ Il998l) . Given their excitation requirements, the 
vibrationally excited SiO maser lines must arise from the 
closest vicinity of the stars. This indicates that high velocity 
components are already present on the smallest scales and 
receive their peculiar high velocities, sometimes exceeding 
the terminal velocity, by other processes than the "regular" 
outflow, which only can commence further out, once dust 
starts forming. A possible mechanism may be related to 
shocks driven into the atmosphere and inner envelope by 
the ar rival of giant conve ction cell s to the star's surfac e 
layers dSchwarzschildll 197*5 ) and, e.g., Jossehn & Plez ( 2007 ). 
Puzzlingly, we see no sign of these extreme velocities in 
any of the water lines, neither in the ground-state lines nor 
in vibrationally excited H2O lines, which all cover much 
smaller velocity ranges than high velocity SiO maser fines, 
although the excitation requireme nts of vibrationally excited 
H?0 and SiO fines a r e similar ([Menten & Melnick Il989t 
Menten & Young! 1 19951 iMenten et alJl2006l) . These lines do, 
in fact, cover velocity ranges comparable to those of the 
321 and 437 GHz H2O lines observed by us. Indeed, under 
plausible physical conditions for the regions from which the 
vibrationally lines arise, n(H2) = 2 10' cm"-' and T = 1000 K 
(iMenten et aL.2006.) . one also finds inversion of these ground 
state fines (see 94.11 1. 

A spectrum of the thermally excited 1 10 - 1 01 ortho ground- 
state H2O line measured with Submillimeter Wave A stronomy 
Satellite (SWAS) has been pubfished by iHarwi t & Bergin| 
(l2002l) . They find a FWZP value of 50 km s similar to that 
we find for the wider of the maser lines. It is, in contrast to the 
case of W Hya, significantly smaller than that of CO lines (see 
above). 

The variety of H2O line shapes toward VY CMa ranges 
from that of the 22.2 GHz spectrum, which entirely consists 
of narrow spikes, to spectra in which most or an appreciable 
portion of the line flux arises from a broad (up to 50 km s"' 
wide) pedestal feature on which a few spikes are superposed 
(in the 325, 439, 471, and 475 GHz fine). The 321 and 437 
GHz spectra are dominated by a few, strong narrow features 
with a little pronounced broad pedestal. The former covers a 
much narrower range than the other lines, just ± a few km s ' 
around the stellar velocity. 

That the kinematics around VY CMa are indeed complex 
also on scales much larger than the masing regions is 
demonstrated b y interferorn e tric C O observations which are 
interpreted by iMuller et al.l (l2007b in terms of spherically 
symmetric outflow together with a high velocity bipolar flow. 



4. Discussion 

4.1. Modeling the excitation of water maser lines 

Earlier excitation studies aimed at explaining stellar (and 
interstellar) H2O maser pumping only considered a limited 
number of energy levels, both, due to limitations of 
computational capabilities and because collision rates were not 
know n for excitation to and from h i gher energy states (see , 
e.g., Ide Jong|[l97l IPeguchil [l977t Icooke & Efitzuil [l985l) . 
Nevertheless, all these studies predicted strong inversion of the 
22.2 GHz line, which was the only H2O maser line known 
at the time. These studies also predicted maser action in a 
number of other transitions that were later indeed found to 
be masing, despite the fact that in some cases the collisional 
excitation rate coefficients use d were vastly different from the 
"modern" values computed bv iPakna et al.l (1988 ). The more 



recent, comprehensive, study by NM91 considered all 349 
states of ortho and para water up to an energy of 7700 K above 
the ground state and extrapolated the Palma et al. rates to high 
en ough temperatures. 

'Neufel d & Mehiickl(ll990l) and NM91 performed statistical 
equilibrium calculations modeling vibrational ground state 
H2O excitation and radiative transfer for plane parallel 
media un der the la r ge vel ocity gradient (LVG) assumption. 
.Neufeld & Mehiickl il990h specifically explore conditions 
under which, both, the the 22.2 GHz and the 321 GHz lines are 
masing for the case of an O-rich evolved star with a mass-loss 
rate (3 10"^ Mq yr ') intermediate between the values of our 
objects. They find that the pumping for these lines is dominated 
by collisions and that radiation plays a minor role. The same 
was also found for other transitions in circumstellar envelopes 
(lDeguchilll977tlCooke & Elitzu^l 19851 NM91). 

NM91 explored larger regions of (essentially) density 
and temperature parameter space in which the energy levels 
of certain transitions become inverted. They calculate maser 
emissivities and optical depths as functions of temperature 
and a variable, which is equal to a geometric factor, G, 
divided by the velocity gradient, dv/dr, times the product of 
the hydrogen nuclei, n, and the water density, n(H20), i.e., 
« X «(H20) = n^Jc(H20), where ;>c(H20) is the H2O abundance. 

is, thus, given by: 



dv/dr 



Assuming a plausible, constant, value for the H2O abundance 
and the velocity gradient of 10""^ and 10"^ cm s"' cm"', 
respectively, makes ^' dependent on the H density squared 
alone. Measuring the H density in units of lO' cm"^ in the 
parametrization of NM91, values of logio(' - -1,0, and 1 
correspond to H densities of 0.32, 1, and 3.2 x lO' cm"^. 

4.2. Comparison with the APEX observations 

NM91 predicted maser emission in all of the lines listed in 
Table ffl except for the 437 GHz 753 - 650 and 443 GHz 
752 - 661 transitions in, roughly, this range of values quoted 
above. Maser emission in all their predicted maser lines has 
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indeed been observed in star-forming regions and/or red giant 
stars, except for the 355 GHz 174 13 - 167,10 transition. They 
present opacities and emissivities (which are proportional to 
isotropic photon luminosities) for assorted lines calculated for 
two values of the kinetic temperature, 400 K and 1000 K. In 
particular, they find the very widespread 22.2 GHz 616 -523 
transition to be by far the most luminous and the most robust 
maser line in the sense that it is inverted over the widest range 
of physical conditions in particular at extreme values of the 
density. A general trend is that for 1000 K the maximum 
emissivities in all the lines are reached for a factor of a few 
higher densities than at 400 K. 

In the calculations of NM91, the peak emissivities of all 
lines but the 22.2 GHz line are similar within a factor of a few 
(although they occur at different values of ('). At, both, 400 
and 1000 K, the peak emissivity of the 22.2 GHz line is much 
higher than that of any of the others, by a factor of ~ 10 at 400 
K and ~ 5 at 1000 K. The peak emissivity is also reached at 
an order of magnitude higher value of than for other lines. 
As can be seen from Table [T] and Fig. |4]in none of our stars 
does the 22.2 GHz line show the behavior predicted by these 
calculations, i.e. that it should be much more luminous than 
all the other lines. While it is among the strongest lines in 
VY CMa, it is anomalously weak in R Leo. We note that R 
Leo showed a n even more extreme 321 /22.2 GHz line ratio 
in 1989/1990 jMenten & Melnicklll99ll) . In this star, the 22 
GHz line is known to undergo extreme and erratic variations 
changing from not or barely detec table to hundreds of Jy within 



dozens of days (iRudnitskijII 1 9871) 



4.3. The "new" 475 GHz transition 

For the newly discovered 475 GHz 533 - 440 transition NM91 
predict, as for the 439 and 471 GHz lines, maser action only 
for their 1000 K calculation. The presence of maser action in 
this line is, thus, an important indicator of a high temperature 
environment. 



4.4. Non- and mispredictions 

NM9ldo not predict maser emission for the 437 GHz 753 - 650 
transition, which we find to be by far the strongest line in 
R Leo and the second strongest in VY CMa. This line has 
only been found tow ard evolved stars and not toward SFRs 
dMelnick et al.lll993h . Pumping by the strong IR field in these 
objects, possibly via line overlaps, not considered by NM91, 
may be responsible for its excitation. NM91 do, however, 
predict maser action in the extremely high excitation (5764.3 
K) 174,13 - 167,10 transition near 355 GHz high with peak 
emissivities at ^' values 1-2 order of magnitude higher than 
values for which other lines show their peak emissivities. For 
constant water abundance and velocity gradient this translates 
into H2 densities of order 10'° " cm"^, which is 1-2 orders of 
magnitudes higher than the values over which most of the other 
lines show maser action. Our non-detection of this line (see 
Table |2] and 34.21) seems to indicates that under such extreme 



conditions the necessary requirements for producing detectable 
maser emission may not be met. 

4.5. A note of caution 

Guided by the described model calculations, observations of 
combinations of maser lines in principle should put constraints 
on the physical parameters of the masing regions. One naive 
hope might be: The more lines one finds masing the tighter the 
physical parameters of the masing region can be constrained. 
One might even hope that the observed relative luminosities of 
the lines might provide yet tighter constraints. 

Unfortunately at present, in the absence of more realistic 
models of the maser regions such expectations are most 
likely overly optimistic for several reasons: Comparing the 
emissivities or their ratios predicted by the generic calculations 
of NM91 for a plane-pamllel(!) medium with measured data 
probably does not make much sense. Moreover, an inherent 
assumption in the values calculated by NM91 is that all the 
lines are saturated throughout the region in which they are 
inverted. There is good observational evidenc e that this is not 
true for the 22.2 GHz line around W Hya ( Reid & MentenI 
ri990) although this issue is unclear (and difficult to address) 
for other transitions and other stars. 

Another caveat when using line ratios as diagnostics is 
time variability. While all our APEX and Effelsberg data were 
taken within 1-2 weeks, at least the 22.2 GHz line in some 
sources (including R Leo) is known to show week-to-week 
variability (see 34.2b . Information on variability on such small 
time scales is not available for the submillimeter lines, but 
presumably all of them are variable. This has been directly 
shown from a comparison of line profiles taken at different 
epoc hs 1-2 months apart for the 321 and 325 (and 22 .2) GHz 
Hnes dMenten & Melnicklll99ll:lYates & Coh"enlll996h 

Furthermore, to use multi-line data as diagnostics, an 
inherent assumption is that all lines arise from the same region 
(whose physical parameters are to be constrained). However, 
for VY CMa, the different shapes of the lines and the different 
velocity ranges they cover suggest that this is not true and 
that at least a significant portion of the lines' emission arises 
from disjoint regions. Clearly, high resolution interferometric 
observations are necessary to get a clearer picture of this. Given 
the high expected brightness temperatures, such observations 
with the Atacama Large Millimeter Array (ALMA) will afford 
studies of R Leo and W Hya with a resolution better than 
a stellar radius a nd a few stellar radii for VY CMa (see 
§4.4 of Menten 2000l) . All the lines discussed here have 
frequencies within ALMA's 275-370 and 385-500 GHz "Firs t 
Light" Receiver bands ("Bands 7 and 8", I Wilson et aDl2005h . 
Observations of the 321 and 325 GHz l ines can alread y be 
made today with the Submillimeter Array jHo et al ] |2004 . 

It is clear that current models for H2O maser excitation 
in circumstellar envelopes fail to explain the observational 
picture. Velocity-resolved observations of a large number of 
maser and non-maser lines from a wide range of energies above 
the ground state with the Heterodyne Instrument for the Far 
Infrared (HIFI) aboard the Herschel satelUte to be launched 
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in the near future will soon deliver very tight constraints on 
H2O excitation. Together with high spatial resolution data for 
high excitation (mostly maser) lines (see above), which will 
only be attainable with ground-based interferometers and more 
sophisticated radiative transfer modeling these observations 
will lead to an understanding of water, which dominates the 
thermal balance in these regions. 

4.6. Different natures of AGB star and red supergiant 
water masers? 



Let us come back to the fact, mentioned in 93.21 that toward 
W Hya, the 22.2 GHz H2O maser emission arises from a ring 
with a diameter of 24 AU centered on the star. In contrast, 
toward VY CMa, the maser emission arises from a 770 x 440 
AU region and proper motion observations ap pear to indicate 



that t he masers partake in a bipolar outflow jRichards et al 



199 A although we notice that this flow would be directed 
perpendicula r to t he larger-scale CO outflow imaged by 
Muller et al.l (l2007h . Are these distributions consistent with 
the masers arising from a "normal", i.e. steadily expanding 
circumstellar outflow? 

In the following, we investigate whether over the whole 
of the maser distributions in both objects the temperatures 
(calculated from the luminosities of the stars) and the densities 
(derived from the mass-loss rates) are conducive for maser 
excitation. 



For W Hya lJusttanont et al.1 (120051) modeled CO data to 
derive a temperature, T, profile, i.e., T vs. distance from the 
star, r. For the maser region (r = 1.8 10'"* cm), one finds 
T ^ 1000 K from their Fig. 3. Also, assuming, like these 
authors, a luminosity of 5400 Lq and an effective temperature 
of 2500 K, respectively, one calculates a comparable 970 K at 
this r from the Stefan-Boltzmann law. 

The H2 density in a spherically symmetric outflow at 
distance r from the central star is given by 

n(H2) = 1.13 10'*^M(Moyr-i)r-2(cm)v;i (kms-') 



where M and Vexp are the mas s-loss rate 
and the expansion velocit y, respectively, (see, e.g. Millar! 1988 : 
Menten & Alcolea 1995 1). If we assume for the masing region 



«(H2) = 10^ cm"^ (see 94.2l l. r - 1.6 lO'"^ cm and an expansion 
velocity of 2 km s"' and invert this equation to calculate the 
mass-loss rate, we obtain M = 6 10"^ Mq yr"\ which is half 
the rate derived by Zubko & Elitzur (2000) an d in betw e en the 



estimates of lNeufeld et all d 19961) and JJusttanont et al.l (120051) 
and, thus, hopefully, a plausible value. 



^ iNeufeld et all h996l) derive a large value of 0.3-2 10"^ Mq yr"' 
for M from modeling emission of thermally excited H2O lines 
observed by the In frared Space Observatory (ISO). In contrast, 
llusttanont et al.l ( l2005i) . also using data from ISO as well as the Odin 
satellite, derive a much smaller mass loss rate of (2.5 + 0.5) lO"' 
Mq yr"', comparable to values derived from CO lines. Their low 
M comes "at the expense" of an extremely high H2O abundance 
of [H2O/H2 = 2 10^, which is significantly higher than the cosmic 
abundance of O. To explain it they invoke an influx of water from 
evaporating icy bodies in an Oort cloud analog, a mechanism that 



Doing the same exercise for the RSG VY CMa with a 
luminosity of ^ 2 10^ Lq, we calculate for a region with 
an average radius of 300 AU (see 93. 2t a temperature of 
480 K, which would be adequate at least for 22.2 GHz 
maser production. However, if we calculate the mass-loss rate 
required for a spherically symmetric outflow to maintain a 
density of 10^ cm"-' at 300 AU, we obtain a value of 0.045 
Mq yr"'. This is « 200 times higher than the generally quoted 
(already extremely high) mass-loss rate for this object (see 
We conclude that, in contrast to W Hya, at least the 22.2 GHz 
H2O masers far from VY CMa cannot arise from the general 
material in a spherically symmetric outflow, but rather from 
special regions (e.g., shock fronts) within such a flow or from 
a coUimated flow. This situation is reminiscent of H2O masers 
in star-forming regions modeled bv lElitzur et al. (Il989h . These 
authors invoke a temperature (400 K) similar to the value we 
derive for VY CMa. 

At this temperature, the highest excitation lines' maser 
emissivities are small (NM91), restricting their emission to 
regions closer to the star. This is in fact consistent with the 
smaller velocity ranges covered by the 321 and the 437 GHz 
lines. 

5. Summary 

We have used the APEX telescope for a multi-transition study 
of the AGB stars R Leo and W Hya and the red supergiant 
VY CMa. We find a variety of relative line ratios between 
the stars. For VY CMa we observe significant line shape 
differences between individual transitions, indicating that the 
lines' dominant contributions arise from different emission 
regions. While for all but one of the maser lines model 
calculation indeed predict inversion in a hot (T = 1000 K), 
dense (n > 10^ cm"^) medium, other predictions of these 
calculations are at variance with the observations. This is not 
surprising, as key assumptions of the calculations do most 
likely not apply. We detected a new, high excitation maser 
line near 475 GHz that is predicted to only arise from a 
hot environment as described above. Studying the excitation 
requirements and the extent of maser emission for W Hya and 
VY CMa, we conclude that (at least the 22.2 GHz) masers in 
the former may be excited in the regular circumstellar outflow, 
while a special excitation mechanism is required for the much 
more widespread emission around VY CMa, possibly shocks 
in a higher velocity outflow. 
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